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ABSTRACT 

Silicates are an important component of interstellar dust and the structure of these grains - amorphous versus 
crystalline - is sensitive to the local physical conditions. We have studied the infrared spectra of a sample of 
ultra-luminous infrared galaxies. Here, we report the discovery of weak, narrow absorption features at 11, 16, 
19, 23, and 28 /im, characteristic of crystalline silicates, superimposed on the broad absorption bands at 10 and 
18 fjm due to amorphous silicates in a subset of this sample. These features betray the presence of forsterite 
(Mg2Si04), the magnesium-rich end member of the olivines. Previously, crystalline silicates have only been 
observed in circumstellar environments. The derived fraction of forsterite to amorphous silicates is typically 
0.1 in these ULIRGs. This is much larger than the upper limit for this ratio in the interstellar medium of the 
Milky Way, 0.01. These results suggest that the timescale for injection of crystalline silicates into the ISM is 
short in a merger-driven starburst environment (e.g., as compared to the total time to dissipate the gas), pointing 
towards massive stars as a prominent source of crystalline silicates. Furthermore, amorphization due to cosmic 
rays, which is thought to be of prime importance for the local ISM, lags in vigorous starburst environments. 

Subject headings: galaxies: ISM — ISM: evolution — infrared: galaxies - infrared: ISM 
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1. INTRODUCTION 

Silicates are an important component of interstellar and cir- 
cumstellar dust. Many Galactic sources show broad emis- 
sion or absorption features at 10 and 18 f im due to the Si-O 
stretching and O-Si-O bending vibrations of silicate bonds. 
Because of the width of these features, these silicates must 
have an amorphous structure. These broad amorphous silicate 
features have also been observed in a variety of extragalactic 
environments, including Seyfert gal axies (e.g. [Laur ent et alJ 
1200(1 IClavel et alJ 12000). quasars llHao et all 120051) lumi- 
nous and ultra- luminous infrared ga laxies ( Genzel et all 19981: 
i RigQDOulou et alJl!999t iTran et alJl200ll Ispoon et aljbooi 
Armus et al. 2004) and a sample of mid-infrared detected, op- 
tically invisible, high -luminos ity galaxies with redshifts of 1 .7 
< z < 2.8 l lHouck et all2005l) . 

In recent years, observational evidence for the presence 
of crystalline silicates in various astrophysical environments 
has also emerged. In particular, infrared spectra have 
revealed that silicates in circumstellar environments often 
contain a significant crystalline fraction around both pre- 
main-sequence stars (T-Tauri stars and Herbig AeBe stars) 
and post-main-sequence stars (Asymptotic Giant Branch 
(AGB) stars, post-AGB stars, planetary nebulae, Lumi- 
nous Blue Variables (LBVs), Red Super Giants (RSGs) 
and post-RSGsl (e g Ia eke et alJl2004t Malfait et1d1ll998t 
[ Molster et all 120021 iSvlvester et all I1999L Ivan Boekel et alJ 
120041 iVoors et all 119991 12000) . These crystalline silicates 
are invariably magnesium rich (e.g., pyroxene (MgSiOs) and 
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forsterite (Mg2SiO/t). Crystalline silicates are also known to 
be ubiquitous in the solar system, including primitive o bjects 
such as comets (e.g.lCrovisier et all 1 9971 lHanner et all 19991: 
iWooden et alJ fT993). Because crystallization is inhibited by 
high-energy barriers, the origin and evolution of the crys- 
talline silicate fraction in interstellar and circumstellar media 
has the potential to provide direct evidence of energetic pro- 
cessing of grains. 

Silicates play an especially large role in shaping the mid- 
infrared spectral appearance of Ultra-Luminous InfraRed 
Galaxies (ULIRGs; L ffi >1O 12 L ). ULIRGs are thought to 
represent the final stage in the merging process of gas-rich 
spirals, where the interaction has driven gas and dust towards 
the remnant nucleus, fueling a massive starburst and a nascent 
active galactic nucleus (AGN). Here, we present evidence for 
the presence of crystalline silicates as part of the deep sili- 
cate absorption features observed towards a sample of twelve 
heavily obscured ULIRGs. 

2. OBSERVATIONS AND DATA REDUCTION 

We are obtaining mid-infrared spectroscopy for a sample 
of 110 ULIRGs, as part of the Guaranteed Time Observa- 
tion ( GTO) ULIRG prog ram of the Infrared Spectrograph 
(TRS) y (Houck et al. 2004) on the Spitzer Space Telescope 
Wern er et all 12004). Seventy seven of these spectra have 
been analyzed so far. 

The twelve ULIRG spectra presented in this paper were se- 
lected based upon the strength of spectral structure indicative 
of crystalline silicate absorption bands at 16 and 23 /im. Ta- 
ble^lists the basic properties of these targets, along with their 
observation dates and on-source integration times. 

The observations were made with the Short-Low (SL) and 
Long-Low (LL) modules of the IRS. The spectra were ex- 
tracted from the flatfielded images provided by the Spitzer 
Science Center (pipeline version SI 1.0.2). The images were 
background-subtracted by differencing the two SL apertures 

9 The IRS was a collaborative venture between Cornell University and Ball 
Aerospace Corporation funded by NASA through the Jet Propulsion Labora- 
tory and the Ames Research Center 
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TABLE 1 
Properties of Sources 



Target AOR key Date observed Int. time redshift Dl 
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00397-1312 


4963584 


04 Jan. 2004 


14 


0.262 


1320 


01199-2307 


4964864 


18 Jul. 2004 


14 


0.156 


737 


06301-7934 


4970240 


1 1 Aug. 2004 


12 


0.156 


700 


06361-6217 


4970496 


1 1 Aug. 2004 


14 


0.160 


760 


08572+3915 


4972032 


15 Apr. 2004 


6 


0.0584 


258 


15250+3609 


4983040 


4 Mar. 2004 


7 


0.0554 


244 


Arp220 


4983808 


29 Feb. 2004 


5 


0.0181 


78 


17068+4027 


4986112 


16 Apr. 2004 


14 


0.179 


858 


18443+7433 


4987904 


5 Mar. 2004 


14 


0.135 


627 


20551^1250 


4990208 


14 May 2004 


4 


0.0427 


186 


23129+2548 


4991488 


17 Dec. 2003 


23 


0.179 


858 



"assuming H =71 kms~' Mpc~' , f2 M =0.27, C A =0.73, fl K =0 

and for LL, by differencing the two nod positions. Spec- 
tra were then extracted and calibrated using the IRS standard 
star HR6348 for S L and the stars H R6348, HD 166780, and 
HD 1735 1 1 for LL (ISloan et al.l2005l) . Small wavelength cor- 
rections were made to compensate for known offsets in the 
S 1 1 processed data. After extraction the orders were stitched 
to LL order 1, requiring order-to-order scaling adjustments 
of typically 5-10%. The largest adjustment was made for 
Arp 220 to match SL order 1 and LL order 2, which required 
SL1 to be scaled up by 21%. In the final step, the 5-37 /jm 
spectra were scaled to match the observed IRS blue or red 
peak-up flux. For those sources lacking a (useful) peak-up 
flux, the spectra were scaled down by 10%, the average scal- 
ing factor of the other spectra. In some spectra, most notably 
in IRAS 20551^+250, residual fringing, which is due to sub- 
pixel pointing errors, appears in LL order 1 (between 20 and 
30 /im). 

3. ANALYSIS 

The top panel of Figure^] shows the average IRS low- 
resolution spectrum for our sample, obtained after scaling the 
spectra to S„ = 1 Jy at 15 /im. The average spectrum is dom- 
inated by broad absorption features at 10 and 18 /im, which 
we attribute to absorption by amorphous silicates. The spec- 
trum further shows weaker absorption features due to wa- 
ter ice (6.0 /im) and hydrocarbons (6.90 and 7.25 /im). PAH 
emission features can be seen at low contrast at 6.2, 7.7, and 
1 1 .2 /im. The features are especially wea k in the spectrum 
of IRAS F00183-7111 ( Sooo netlflEOOl and absent in the 
spectrum of IRAS 08572+3915 (Figure^ lower panel). Our 
spectra further show a few emission lines, most notably the 
H 2 S(3) line at 9.66 /im, the [Nell] line at 12.81 /zm and the 
FL- S(l) line at 17.0 /zm. These will b e discussed elsewhere 
( Armu s et al.l2005l Iffi gdon et al.l2005l) . 

Upon close inspection, the spectra also reveal evidence for 
weak and narrow absorption features near 16, 19 and 23 /zm, 
independent of the redshift of the source. Emission and ab- 
sorption features at these wavelengths are also know in Galac- 
tic sources with strong silicate features and have been at- 
tributed to the presence of crystalline silicates. 

In order to investigate the presence of a crystalline com- 
ponent to the silicate absorption features in our sample, we 
infer the silicate optical depth spectrum by dividing out a 
local continuum from our IRS spectra. We define the lo- 
cal continuum as a three point spline interpolation of contin- 
uum points at 5.6 /zm (0.1 /zm shortward of the wavelength 
range affected by the 6.0 /zm water ice absorption feature), 
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FIG. 1 . — Upper panel: the average 5-32 (im spectrum of the 12 ULIRGs 
in our sample after normalizing each spectrum individually at 15 /im. Fea- 
tures and lines are identified in the text. Lower panel: the 5-60 fim spectrum 
of IRAS 08572+3915. The IRS low-resolution spectrum (black) is shown to- 
gether with the 12, 25 and 60 /im IRAS fluxes (grey filled circles). Three 
choices for the local continuum are indicated by the black dashed line and 
grey shaded area. Error bars are ommited in both panels, because of the high 
S/N of the data 



7.1 /zm (in between the two hydrocarbon absorption bands 
at 6.90 and 7.25 /im), and the red cut-off of IRS LL order 
1 (ranging from 29 to 36 /zm, depending on the redshift of 
the source). For sources without discernable 7.7 /zm PAH 
emission (e.g. IRAS 08572+3915 and IRAS F00183-7111), 
we replace the 7.1 /im pivot by the continuum at 7.9-8.0 /zm. 
The resulting spline interpolated local continuum is illustrated 
for IRAS 08572+3915 in the lower panel of Figure^ The 
effect of small changes in the adopted local continuum on 
the resulting silicate profile is illustrated by the shaded areas 
around the spline-interpolated continuum in the lower panel 
of Figure^ and on the resulting silicate profile in the third 
panel of Figure|2] The resulting optical depth profiles for 
the twelve sources in our sample have apparent 10 /zm opti- 
cal depths ranging from 2.1 for IRAS 17068+4027 to 4.2 for 
IRAS 08572+3915, with a mean value of 2.5 (see Figure|3]and 
Table|2j. 

In order to derive the characteristics of these narrow fea- 
tures we have to account for the amorphous silicate features 
which dominate the absorption structure in this wavelength 
range. We have compared the general profile of the 10 and 
18 /zm amorphous silicate featues in these ULIRG S with those 
in a s ample of Galactic background sources ( Chia r & Tielensl 
120051) . However, because of the uncertainties in the placement 
of the continuum for these Galactic sources, the derived inter- 
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FIG. 2. — Illustration of our method for fitting the amorphous compo- 
nent of silicate absorption features. Upper panel: laboratory opacity spectra 
of amorphous olivine (black) and amorphous pyroxene (grey), presented as 
optical depth profiles. Overlayed are spline fits (red, blue) to pre-selected 
wavelength pivots (red, blue diamonds). Fit residuals, indicating imperfec- 
tions in the fitting method, are shown as deviations from zero optical depth 
(red, blue). Second panel: spline fit (red) to the GCS 3 optical depth pro- 
file (black; scaled-up by a factor 1.75). Fit residuals are shown as deviations 
from zero optical depth (red). Third panel: the optical depth spectrum of 
IRAS 08572+3915 (black). The grey shading indicates the effect of a dif- 
ferent choice of local continuum on the resulting optical depth profile; see 
FigurefTI The red curve is a spline fit to the amorphous silicate profile, ignor- 
ing the narrow substructure near 11, 16, 19 and 23 /im. The residual spectral 
structure, shown in red at the top of this panel, is too strong to be attributed 
to imperfections in the fitting method. Instead we attribute the features at 1 1 , 
16, 19 and 23 fim to crystalline silicates. Third panel inset: a close-up of the 
14—21 fim range, demonstrating the prescribed placement of spline pivots at 
14.5, 15.0, 17.0, 17.5, 20.0 and 21.0 /an. Lower panel: the residual optical 
depth spectrum of IRAS 08572+3915 (red is compared to the residual crys- 
talline silicate emissivity of the disk around the young star IS Tau (blue) and 
the opacity profile of forsterite (purple). Both have been scaled by arbitrary 
factors 



stellar 10 and 18/im amorphous silicate features show differ- 
ences in their relative strengths, and in the level of absorption 
between them. Since fitting the ULIRG data with a compos- 
ite Galactic spectrum would therefore introduce spurious ar- 
tifacts, we have elected to represent the broad amorphous sil- 
icate features in the ULIRG and Galactic spectra with spline 
fits. Deviations from these fits in the ULIRG spectra, sig- 
nificantly larger than those seen among the Galactic sources, 
would then argue strongly for a crystalline component. We 
have tested this procedure on the high signal-to-noise spec- 
trum of the Galactic Center Source GCS 3, and on laboratory 
absorption spectra of amorphous silicates. The residuals from 
the spline fits to these spectra are very small (<6%), validat- 
ing our analysis to that level. 

The method is demonstrated in the top panel of Figure|2] 
for laboratory spec tra of amorphous olivines and pyroxenes 
(Fabian et all l2001l) . The spline fit traces the amorphous pro- 
file to better than 2% of the local optical depth, except for 
the 16 /jm range, where the spline deviates by up to 6%. 
The resulting small artifacts are shown at the top of the 
panel. We also tested our method on the optical depth pro- 
file as seen towards th e Galactic background source GCS 3 
( Chiar & Tielensll2005h . depicted in the second panel of Fig- 
ure|2] The silicates along this line of sight are thought to 
be > 99% amorphous in composition ( Kemper et alJ 120041 
2005). In contrast to the laboratory profiles, the fit residu- 
als for the GCS 3 spectrum are dominated by spectral noise 
rather than fitting artifacts. The third panel of Figure[2]shows 
the fit to the silicate profile of the ULIRG IRAS 08572+3915. 
The red curve smoothly fits the maximum depth of the broad 
18 /im feature, while ignoring the strong, narrow substructure 
around 11, 16, 19 and 23 /im. After subtracting the spline- 
interpolated amorphous component from the optical depth 
spectrum, the features at 11, 16, 19 and 23 /im show up in the 
fit residual spectrum, clearly above the levels expected for ar- 
tifacts introduced by the fitting procedure. In the lower panel 
of Figure[2]we compare these residuals to known spectra of 
circumstellar and laboratory crystalline silicates. 

The adopted spline fits to the amorphous silicate profiles 
for the twelve ULIRGs in our sample are overplotted in red in 
Figure|3] while the residual optical depth spectra are presented 
in Figure^] The individual spectra in the latter figure are trun- 
cated at -0. 1 optical depth in order not to be dominated by spu- 
rious optical depth structure introduced by emission features 
of 1 1.3 /im PAH and 12.8 /im [Ne II]. These common features 
aside, the twelve spectra clearly show absorption structure at 
16, 19 and 23 /im that cannot be attributed to observational ar- 
tifacts (i.e. fringing, flux calibration) or imperfections in the 
spline fit to the amorphous component. A few sources (e.g. 
IRAS 06301-7934, IRAS 23129+2548) further show a weak 
feature at 27-28 /im and only one source, IRAS 08572+3915, 
a strong feature at 1 1 /im (see also Figure|2j. The latter fea- 
ture may also be present in the other spectra, but its detection 
is complicated by the presence of emission from 10.51 /im [S 
IV] and 11.3 /im PAH in the same wavelength range (see Fig- 
ure0}. Further note the presence of absorption bands of gas 
phase C2H2 and HCN at 13.7 and 14.05 /im in the spectra of 
several of our ULIRGs, most notably IRAS 15250+3609 and 
IRAS 23129+2548. These bands are otherwise seen in dense, 
warm molecular clouds surrounding massive Galactic proto- 
stars (Lahuis & van Dishoeck 2000). 

In Figure|5] we show the average residual optical depth 
spectrum for our twelve sources and compare it to the ob- 
served crystalline silicate opacity profile of the young star IS 
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FIG. 3. — 8-30 fim optical depth spectra for the 12 ULIRGs in our sample. The red curve is a spline fit to the amorphous silicate profile, ignoring the narrow 
substructure at 11, 16, 19. 23 and 28 fim. Vertical dotted lines indicate the positions of 9.66 fim H? S(3), 10.5 /im [SlV], 1 1.3 fim PAH, 12.8 fim [Nell], 15.5 /im 
[Nem], 17.0 fim H 2 S(l) and 18.7 /im [Sill]. 



Tau JSargent e t al.l200 5|) and the calculated opacity profile of 
forsterite ( Fabia n et al.l2001l) . The relative strengths of the ab- 
sorption bands in the ULIRG spectra are different from those 
of the comparison profiles. In the ULIRG spectra, the features 
get progressively weaker towards longer wavelength, whereas 
no such trend exists for the comparison profiles (see also Fig- 
ures |21 and [5}- A possible explanation for this trend could be 
dilution of the crystalline absorption spectrum by cold dust 
emission from a less obscured cooler component, either the 
other nucleus or the circumnuclear environment. Neverthe- 
less, the similarities between the ULIRG profile and the com- 
parison profiles in peak position and width are striking. We 
therefore conclude that the features observed at 11, 16, 19, 



23 and 28 /im in the ULIRG spectra are indeed caused by the 
presence of crystalline silicates in their nuclear medium. 

The peak position of the 16 ba nd, in particular, is sensi- 
tive to the Mg/Fe ratio of olivines (Koike et al. 2003). The 
peak position of this band observed in the ULIRG spectra 
(16.1 /im) falls close to the laboratory measured position of 
forsterite . the Mg-ric h end member of the olivines (Mg2Si04; 
16.3 um: lKoike et al.l2003l) and much to the blue of the peak 
position in fayalite, the iron-rich end m ember of the olivines 
(Fe2Si04; 17.7 fim; Koike et all 120031) . Hence, in line with 
other observations of crystalline silicates in space, extragalac- 
tic olivines appear to be extremely Mg-rich and Fe-poor. 

We infer the fraction of crystalline silicates in our sample 



5 




Rest wavelength [/im] 

FIG. 4. — 10-30 fim residual optical depth spectra for the 12 ULIRGs in our sample after subtraction of the amorphous silicate component from the optical 
depth spectra. For plotting purposes, the spectra have been offset and truncated at residual optical depths of -0.1. Crystalline silicate features can be identified at 
11, 16, 19, 23 and 28 /xm 



TABLE 2 

Amorphous and crystalline silicates 



Target 


T fl „,(peak) a 


T„-(peak) b 


r cr (peak) lr m 


(peak) N„-/N,-„„ c 


00183-7111 


3.0 


0.22 


0.073 


0.11 


00397-1312 


3.2 


0.28 


0.088 


0.13 


01199-2307 


3.6 


0.21 


0.060 


0.09 


06301-7934 


3.4 


0.35 


0.10 


0.15 


06361-6217 


2.9 


0.15 


0.051 


0.08 


08572+3915 


4.2 


0.22 


0.053 


0.08 


15250+3609 


3.6 


0.20 


0.056 


0.08 


Arp220 


3.3 


0.16 


0.048 


0.07 


17068+4027 


2.1 


0.17 


0.082 


0.12 


18443+7433 


3.2 


0.25 


0.078 


0.12 


20551^1250 


2.9 


0.18 


0.062 


0.09 


23129+2548 


3.2 


0.25 


0.078 


0.12 



"apparent peak optical depth of the 10 (im feature 
b apparent peak residual optical depth of the 16 fim feature 
c assuming Ncr/Nnm = 1.5 X T cr (peak) /r am (peak) ; see eq. (1) 

from the peak optical depths of the 10 /im amorphous and the 
16 /im crystalline silicate absorption bands using 



Ncr_ 

Nam 



r cr (peak) /c flffl (peak) 

X 



(peak) K cr (peak) 



(1) 



where N am and N cr are the mass column densities of the 
amorphous and crystalline silicates. Adopting K I/m (peak) = 
2.4xl0 3 cm 2 /g as th e peak mass absorption coefficient for 
amorphous silicates (iDorschner et alJ Tl995) and K cr (peak) 
= 1.6 xlO 3 c m 2 /g as the peak m ass absorption coefficient 
for forsterite JFabian et alJ 12001 ). we find crystalline-to- 
amorphous silicate mass column density ratios ranging from 
0.07 to 0.15, with a median value of 0.11 (see Tabled). It 
should be understood that these numbers are upper limits, 
since foreground emission and radiative transfer effects within 
the optically thick 10 /im silicate band may cause the apparent 



optical depth of the 10 /im silicate feature to be a lower limit 
to the true optical depth. Our findings should be contrasted 
to the upper limit of cryst alline silicates in the gene ral ISM 
of the Milky Way of < 1 % ( Kemper et a lJ2004tl2005l) and the 
crystalline-to-amorphous ratios of up to 0.75 observed in cir- 
cumstellar environments (Molster et al. 2002). 

Given the broad range in infered crystalline-to-amorphous 
silicate ratios in our sample, we have investigated the exis- 
tence of correlations with the apparent 10 /im optical depth, 
the IRAS R(60,100) color, the 5.5 /im-to-23 fjm rest frame 
spectral slope and the 6.2 /im PAH equivalent width, but 
found no significant trends. The fraction of crystallinity 
within this sample can therefore, at present, not be linked to 
the infrared spectral appearance. 

4. DISCUSSION 

We have detected crystalline substructure in the silicate ab- 
sorption features towards twelve strongly obscured ULIRG 
nuclei. This is the first detection of crystalline silicates 
in any source outside the Local Group. Space-based and 
ground-based observations have revealed the presence of 
infrared emission features of crystalline silic ates around 
young stellar obje cts and evolved sta rs jAc ke et alJ 12004: 
Malfa it et all 1 9981 iMolster et al.ll2002tlSvlvester et al.lll999t 
van Boekel et all2004HVoors et al.ll999t 120001) . Comets also 
ofte n show evidence for crystalline silicates in their spec- 
tra (|Crovisier et all 1 19971 lHanner et alJll999t IWooden et alJ 
1993). In contrast, while stellar sources of silicate dust inject 
at least 5% of their silicates in crystalline form, the Galac- 
tic interstellar silicate feature is exceedingly smooth, and 
there is no evidence for a crystalline absorption component 
in the Galactic ISM. This translates into an upper limit on the 
crystalline fraction in the Galactic interstellar medium of 1% 
( Kem per et aDl2004l l20Q5h . This difference in crystallinity 
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FIG. 5. — The average 13-30 fim residual optical depth spectrum for our 
sample (black) compared to the residual crystalline silicate emissivity of the 
disk around the young star IS Tau (blue) and the opacity profile of forsterite 
(purple). Both have been scaled by arbitrary factors 

between silicates injected and silicates in the ISM implies a 
rapid transformation of crystalline silicates into amorphous 
silica t es in t he interstellar medium (< 10 8 yr; Kemper et al. 
2004 l200l . This transformation has been attributed to 
energe tic processin g of the dust by heavy cosmic ray ions 
iBringa et al.l 120051) or to i on bombardment i n high veloc- 
ity (y > 1 000 km/s) shock s dJager et all 120031 iRotundi et all 
2002; Demvketal. 2004). In particular, based upon labo- 
ratory experiments and estimated cosmic ray fluxes in our 
galaxy, the timesca le for am orphization is estimated to be only 
70 million years (Bringa et all 120051) . considerably shorter 
than the timescale at which (crystallin e) silicates are in jected 
into the Galactic ISM (4 billion years; Brin ga et al.l2005l) 

These same processes — injection of crystalline silicates 
from stars and cosmic ray and shock amorphization — will 
play a role in the interstellar media of these ULIRGs as well. 
Possible explanations for the much higher crystalline silicate 
fraction in these ULIRGs as compared to the Milky Way in- 
clude a higher fraction of crystalline silicates injected and/or 
a delay in the interstellar amorphization rate. 

In particular, ULIRGs are characterized by a high rate of 
star formation driven by merging events. In contrast to the lo- 
cal ISM, the enrichment of the dusty interstellar medium will 
be dominated by the rapidly evolving massive stars and the 
contribution by the more numerous low-mass stars will lag on 
the tim escale associated with the U LIRG-merger event (10 8 - 
10 9 yr; Murph vet al.lll996t 12001). Evidence for the pres- 
ence of a large population of evolved massive sta rs in ULIRG 
nucle i exists for the nearest ULIRG, Arp220 iArmus et al.l 
119951) . Spectroscopic studies of sources in our galaxy have 
revealed that massive stars are prominent sources of crys- 
talline silicates during the supergiant phase. In particular, 
the more extreme examples of this class, such as AFGL4106, 
NMLCyg and IRC + 1 0420, have cry stalline silicate fractions 
of the order of 0.15 l lMolster et al.l lT999). Likewise, crys- 
talline silicates have been observ ed in the ejecta of som e 
LBVs, such as R71 and AGCar l lVoors et al.lfl999l 120001 
This episodic mass-loss phase dominates the dust injection 
by the most massive stars (M>50M Q ) in the Galaxy. Finally, 
it is currently unknown how much crystalline silicate dust 
is injected by supernovae. Indeed it is unknown how much 
dust is injected by SNe. However, it is conceivable that the 
SNe resulting from this starburst also inject a large fraction 
of their freshly condensed dust in the form of crystalline sili- 
cates. Hence, a starburst may rapidly increase the crystalline 



silicate fraction in the nuclear region. 

Of course, the high supernova rate in a starburst will drive 
strong shock waves into the local environment, which will 
sputter the dust, returning the a toms from the solid to the gas 
phase dTones et al]ll994l [jj?96l|). However, likely, the destruc- 
tion rate of crystalline and amorphous silicates are similar 
and hence, to first order, this will not affect the crystalline- 
to-amorphous fraction. More importantly, supernovae are 
thought to be the dominant source of cosmic rays and this may 
increase the rate at which crystalline silicates are transformed 
into amorphous silicates. However, cosmic rays freshly accel- 
erated in the supernova remnants ma y leak away to t he rest of 
the galaxy on a timescale of ^10 6 yr ( Webbe rll993l) . limiting 
the rate of this solid-state conversion in the starburst environ- 
ment. Thus, the crystalline-to-amorphous conversion time- 
scale may actually be quite similar to the value cal culated for 
the local Milky Way (70 million years; Bringa et al. 2005). 
We expect that over time, as the merging event ages and the 
starburst intensity decreases, the fraction of crystalline sili- 
cates will decrease on a similar time-scale. Thus, we attribute 
the high fraction of crystalline silicates in these ULIRGs 
as compared to others to the relative 'youth' of these sys- 
tems; e.g., the amorphitization process may lag the merger- 
triggered, star-formation-driven dust injection process. 

All 17 ULIRGs with r(10^m)>2.9 in our ULIRG sam- 
ple of 77 galaxies show crystalline silicate structure in their 
IRS spectra. At lower silicate optical depth, the number of 
clear detections falls sharply. Above r( 10 /im)=2, 18 out of 23 
ULIRGs definitely show the 16 crystalline feature. Above 
r(10/xm)=l, this is 20 out of 46 and for the entire ULIRG 
sample only 21 out of 77. We note that a complete analysis of 
this fraction is hampered by possible radiative transfer effects 
in the silicate features (as evidenced by the strong variation in 
the ratio of 10-to-18 peak optical depths), the presence of 
16-18 fim UIR emission in some sources with non-negligible 
PAH emission ( Brandl et al. 2005) and by the challenge to de- 
tect weak 16 /im absorption features at low contrast and low 
signal-to-noise. Nevertheless, it is clear that crystalline sil- 
icates are a common component of the nuclear interstellar 
medium during the strongly obscured evolutionary phase of 
the ULIRG phenomena. 

Crystalline silicates can also be a sign of high-temperature 
grain processing. In a ULIRG, both a central AGN and hot, 
young stars are potential sources of high-energy photons ca- 
pable of heating the dust in the ISM. X-ray and far UV pho- 
tons from a central AGN will produce a hot, inner region. 
However, an AGN as the source of grain processing in our 
ULIRG sample is problematic for two reasons. First, interfer- 
ometric 10 fim spectra of the Seyfert-2 nucleus of NGC 1068 
do not show any indicati on for cry stalline silicates in its in- 
ner warm 2pc region (Jaffe et al. 2004). Second, the crys- 
talline silicates are seen in absorption in our spectra, imply- 
ing that they are located in the cooler, outer regions. This 
would require a large scale mechanism to transport the inner 
toroid material outwards and distribute it over the surrounding 
medium. The extended nature of a starburst does not require 
such a mechanism and therefore seems a more likely source 
of processed material than a central source. 

Crystalline silicates are also a characteristic of the cir- 
cumstellar planetary disks surrounding young stellar objects 
such as Herbig AeBe stars and T-Tauri stars. Again, in or- 
der to be seen in absorption, this material will have to be 
transported outwards into the surrounding cold m edium, pre- 
sumably through jets and winds (Tielens 2003). There is, 



7 



however, no indication in Galactic sources for such large 
scale transport of crystalline silicates. On the contrary, for 
Herbig AeBe stars the emission features of crystalline sil- 
icates are stron gly concentrated toward s the inner 2AU of 
the disks (e.g., Ivan B oekel et al. 2004). Hence, we deem 
high-temperature crystallization of existing amorphous sili- 
cate grains unlikely as the source of the high crystalline sili- 
cate fraction in ULIRGs. 

5. CONCLUSIONS 

In this work, we report the discovery of crystalline sub- 
structure at 11, 16, 19, 23 and 28 /mi in the amorphous sil- 
icate features towards twelve deeply obscured ULIRG nuclei. 
These features indicate the presence of the mineral forsterite 
(Mg2Si04). Previously, crystalline silicates have only been 
observed in circumstellar environments. 

We infer the fraction of crystalline silicates in our sample 
from the peak optical depths of the 10 /im amorphous and the 
16 /im crystalline bands and find a crystalline-to-amorphous 
ratio ranging from 0.07 to 0.15, with a median value of 0.11. 
These numbers are likely upper limits, since foreground emis- 
sion and radiative transfer effects will cause the apparent 
10 /im silicate optical depth to be a lower limit to the true 
optical depth. 

The crystalline-to-amorphous ratio in our twelve deeply ob- 
scured ULIRGs is 7-15 times larger than the upper limit for 
this ratio in the interstellar medium of the Milky Way. This 
suggest that the timescale for injection of crystalline silicates 
into the ISM is short in a merger-driven starburst environment 



(eg., as compared to the total time to dissipate the gas), point- 
ing towards evolved massive stars (red supergiants, LBVs 
and type II supernovae) as prominent sources of crystalline 
silicates. Furthermore, the timescale for amorphitization of 
crystalline silicates, which is known to be fairly rapid in the 
ISM of the Milky Way (~ 10 8 yr), is at most of similar or- 
der in starburst environments. We expect that over time, as 
the merging event ages and the starburst decreases in inten- 
sity, the fraction of crystalline silicates will decrease rapidly. 
Thus, we attribute the high fraction of crystalline silicates in 
these ULIRGs, as compared to others, to the relative 'youth' 
of these systems. 

Finally, other galaxy types may also exhibit crystalline sili- 
cate features far above the upper limits set for the crys talline 
silica t e frac tion in the ISM of our galaxy (< 1 %; Kemper et al. 
2004, 2005). The possible detection of a 23 /im crystalline 
emission feature in the quasar PG 1351+640, reported by 
lHao et al 1 d2005l) . may be a signpost of more detections to 
come. 
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